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We propose a phenomenological model where the gravitational interaction between dark matter 
and baryons is suppressed on small, subgalactic scales. We describe the gravitational force by adding 
a Yukawa contribution to the standard Newtonian potential and show that this interaction scheme 
is effectively suggested by the available observations of the inner rotation curves of small mass 
galaxies. Besides helping in interpreting the cuspy profile of dark matter halos observed in N-body 
simulations, this potential regulates the quantity of baryons within halos of different masses. 

PACS numbers: 98.35.+d, 98.62.Gq, 95.30.Sf 



Lacking evidence of direct detection, the presence, na- 
ture and quantity of dark matter (DM) must be in general 
inferred from the kinematic and distribution properties of 
baryons. A model based on the growth of small fluctua- 
tions through gravitational instability in a universe domi- 
nated by cold dark matter (CDM) provides a satisfactory 
fit to a wide range of properties of visible structures on 
scales > 1 h^^Mpc Recently, however, some apparent 
inconsistencies in relating DM and baryon properties on 
subgalactic scales have been highlighted, prompting the 
investigation of radical models to explain the discrepancy. 

Low surface brightness (LSB) galaxies are dominated 
by dark matter ^ whose distribution can be traced by 
the rotation properties of a dynamically unimportant 
stellar component. Therefore, these galaxies provide the 
best laboratories where to test the consistency of the 
CDM scenario. Recent high-resolution studies of the in- 
ner (< 5 kpc) rotation curves of LSB galaxies suggest 
that the dark matter is distributed in spherical halos with 
nearly constant density cores (< 1 kpc) 0,1113. These 
observations seem to collide with the predictions of CDM 
models which produce halos with compact density cusps 
and steep, universal mass-density profiles p ~ r^'^ with 
f3 ~ 1-1.5 (e.g. This disagreement must now be 

taken seriously since it is not easily interpreted in terms 
of observational uncertainties [sl or artificial resolution 
effects in the CDM simulations '9^ . 

A parallel problem is that hierarchical clustering the- 
ories with scale-invariant primordial perturbations show 
significantly more virialized objects of dwarf-galaxy mass 
in a typical galactic halo than are observed around 
the Milky Way Possible solutions to this conflict 

range from non-standard models of inflation |lll | or mod- 
els in which the Universe is dominated by warm, self- 
interacting or annihilating DM (e.g. ^3)' hydrody- 
namical mechanisms that prevent accreting baryons to 
cool into stars (e.g. 0). 

In the absence of an obvious astrophysical mechanism, 
the purpose of this note is to investigate if the various 
mismatches between observations and CDM predictions. 



which however appear to be characterized by a common 
physical length at subgalactic scales, reflect a more fun- 
damental theoretical problem concerning the nature of 
the gravitational attraction between exotic (dark) and 
standard (visible) matter. 

Occasionally, modifications of gravity have been pro- 
posed on scales larger than those investigated here in or- 
der to describe the fiatness of the rotation curves without 
the need of dark matter [3 (see |0| for a review). The 
picture that we propose, on the contrary, aims to mod- 
ify neither the successful DM paradigm nor the standard 
baryonic physics at any scale length. We speculate that, 
in the Newtonian limit of approximation (small curva- 
tures and small velocities in Planck units), the gravita- 
tional interaction in the mixed, dark-visible sector (i.e. 
that between baryonic and non baryonic particles), is 
dampened on small kpc) scales. 

We parameterize a general small-scale modification of 
gravity in the mixed sector, by adding to the standard 
"l/r" potential a Yukawa contribution which is active 
only between visible and dark matter particles. This 
choice is motivated by the known short-range character 
of the Yukawa interaction and by its simple analytical 
properties in both position and momentum space. 

Accordingly, we rewrite the total gravitational poten- 
tial acting on the baryons as 4> = 4>n + ct 4>y , where ip^ 
is the usual Newtonian potential, 

V^^AT = 47:G{pd+Pb), (1) 

and where, weighted by the strength parameter a, 4>y is 
the Yukawa contribution: 

(V2-A-2)0^=47rGpD. (2) 

The parameter A is the typical length scale above which 
(f)Y dies out exponentially. Note that the Yukawa con- 
tribution to the total gravitational potential "felt" by 
baryons is sourced by the DM energy density p^ only. 

By inverting eq. 10) the familiar exponential behav- 
ior of the Yukawa potential is recovered. Neglecting the 
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contribution of the baryons in eq. (0, the gravitational 
potential (?I)(x) felt by a visible test particle in the pres- 
ence of a distribution of DM p£)(x') amounts to 



(x) = -G d'x'pui^r—^^ 
J |x - x' 



(3) 



For the critical value a = — 1 the gravitational attraction 
between two point particles of different species does not 
diverge and is maximally suppressed in the small distance 
limit. We have tested our hypothesis of gravitational sup- 
pression by fitting the parameters of this Yukawa model 
using the rotation curves of LSB galaxies. 

We perform the integration ^ in the case of a spheri- 
cally symmetric distribution of DM and we find that the 
contribution of a spherical shell of radius r' and thickness 
dr' to the Yukawa part of the potential is given by 



d(t>Y{r) 
2ttG 



z ^ " I r pD\r )dr . 



Note that since Gauss theorem does not hold for (see 
eq. (0)) a spherical shell at r' generally does contribute 
to the net force at some internal point r < r' . 

A general prediction of CDM models is that dark mat- 
ter halos should have a steep central cusp, i.e. their den- 
sity profiles have pure power-law forms. By introducing 
a rescaled dimensionless radius x = r/X, we thus focus 
on a density profile of the type Pd{^) = Po . In this 
relevant case the integral Q can be solved analytically in 
terms of confluent hypergeometric functions of the first 
kind $(a; &; z). The velocity curve, = r\d(j)/dr\., is 
given by 



47rGpoA2 
where 
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[F{x,p)+F{~x,l3)] 



(4) 



F{x,P)= 1-- 



(l-^) [N2-/5<i>(l;3-/3;^)-r(3-/3)e 



At large radii [x ^ l,r A), the velocity curve 
Q flattens down to the standard Newtonian behavior, 
V oc a;*-^"^^/^. In the inner regions, however, the Yukawa 
contribution is efficient and the curve may be approxi- 
mated by a different power-law; in this case the Newto- 
nian exponent for the velocity gets a correction that we 
have estimated as —a[\ -\- /3(1 — a)]/ll. In the case of a 
NFW 7] inner profile (/? = 1), eq. Q reduces to 



v'^{x) 
47rGpoA2 



+ -[l-e--(l + x)], (/3=1). 



We have investigated the consistency and the univer- 
sality of this gravitational paradigm, by fixing the model 
parameters by means of the observed (but theoretically 
unexplained) kinematics of low surface brightness galax- 
ies. To this purpose we have used high spatial resolution 
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FIG. 1: Velocity models derived by fitting the high resolution 
rotation curves of LSBC F571-8 |||, NGC 4325 [f, NGC 3109 
fT3|, and NGC 4605 [13 with eq. 5. The fit is performed 
assuming no dynamical contribution from the disk (minimal 
disk hypothesis) . The best fitting parameters for each galaxy 
are shown in the corresponding panel. Errorbars represent la 
uncertainties. 



Hq, rotation curves measured and reduced by different 
authors 0, 0, 0, 0| according to independent observa- 
tional strategies and technics (see fig. 1). 

We have assumed that the LSB halo is the dominant 
mass component and we have not attempted to model 
and remove the contribution from the galactic disk to 
the rotation, which is taken to be negligible (minimal 
disk hypothesis, i.e. the mass-to-light (M/L) ratio of the 
material in the disk is M/L ~ 0, e.g. [l3)- Our results 
may be described as follow: 

a) the parameter a is remarkably stable for all the objects 
considered. Moreover, its mean value a ~ —1 ±0.1 phys- 
ically suggests that, independently of the phenomenolog- 
ically adopted Yukawa model, the effective force between 
baryons and DM tends to vanish on scales smaller than 
some characteristic length A; 

b) the best fitting value for (3 (1.35 ± 0.05), the inner 
slope of the halo mass density universal profile, is in ex- 
cellent agreement with the CDM calculations by 0,0- In 
other terms the observed kinematics of visible matter in 
dwarf galaxies is dynamically induced by the dark mat- 
ter clumps predicted by numerical simulations. As it is 
well known, a pure Newtonian potential fails to reconcile 
observations with theoretical expectations. For example, 
observations would suggest bimodal DM profiles, i.e. de- 
viations from the predicted single power law distribution 
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in the core of the DM clumps (e.g. Even assum- 

ing that there are at least two distinct mass components 
(halo and disk) that are contributing to the Newtonian 
gravitational potential jl7| , a single power law model may 
be recovered but at the cost of a significantly flatter slope 
{P ~ 0.6) than predicted by simulations; 
c ) the scale below which our proposed modification of the 
Newtonian DM-baryon interaction becomes effective is 
less robustly constrained. While the data of and 
seem to point, consistently, to a value A = f .1 ± 0.08 kpc 
(assuming Hq = 70 km s~^ Mpc~^), the rotation curve 
of NGC4605 ([13) favors a smaller value (A - 0.2 ±0.02) 
kpc. Future velocity surveys of the inner regions of lo- 
cal dwarf galaxies will help in better constraining the 
range of variation of this fundamental length and in un- 
derstanding if, as we assume here, NGC 4605 may be 
considered as non representative of the general kinemat- 
ics of dwarf galaxies. Finally we note that above this 
kiloparsec scale, in the external halo regions, a standard 
Newtonian regime is recovered, and a NFW density well 
describe the observed flattening of the rotation curves. 

We have also investigated how this proposed universal 
form for the coupling of the DM and baryons may help in 
interpreting some CDA^I results. Many authors (e.g []J|) 
have argued that in order to match CDM models with ob- 
servations, some mechanism must prevent a large number 
of dwarf galaxies from forming in low mass halos. Much 
attention has been focused on investigating astrophysical 
and hydrodynamical mechanisms that may suppress the 
efficiency of conver ting baryons into stars in small galac- 
tic systems (e.g. !2flj): in particular supernova-driven 
winds that expel a large fraction of the baryonic compo- 
nent from a dwarf galaxy (e.g. 0|, but see HI]), and 
photoionization effects that prevent baryons to cool into 
stars (e.g. 23j ) . 

In our picture, the lack of a significant luminous com- 
ponent in small halos may be interpreted noting that, in 
addiction to local, baryonic physical processes, a cosmo- 
logical mechanism contributes to regulate the quantity 
of baryons that falls into halos of different masses. The 
growing of the baryonic density contrast Sb — Sps/ps 
after decoupling is described in linear perturbation the- 
ory, neglecting radiation effects [23. |25| , by the equation 

Sb + 2HSb = V2(50, (5) 

where a dot means differentiation with respect to the 
proper time i, H ^ d/a is the bubble parameter and S(j) 
is the perturbation of the gravitational potential. DM 
perturbations, approximately the only gravity sources, 
grow almost independently 6d = Spo/po — {a/adec)S']f'^ 
and create the potential wells where baryons fall in. 
In the standard scenario V^dcf) ~ AnGSpo, and the 
RHS term of lO very efficiently drives the growth of 
baryonic perturbations of any length scale (the fraction 
5b /5d — 1 — idec/a, in fact, becomes of order unity 
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FIG. 2: Upper: time evolution of the fraction Sb/Sd pre- 
dicted by our model {a — — 1,A — 1.1 kpc) for perturba- 
tions corresponding to galaxies of different masses. The scale- 
invariant standard result is also plotted (thick line). The 
shaded region represents the post-turn-around epoch. Lower: 
Baryon fraction at the turn-around epoch in units of its mean 
cosmological value Qb /i^o) as a function of the mass of the 
collapsed object and for different values of A. 

already at redshift, say, z — 100). By using Q and 101, 
on the contrary, the driving term on the RHS of (|SJ) be- 
comes scale-dependent: V'^Scj) = 4:TtG{1 — X'^V'^)~^SpD- 
We can gain some insight by noting that, as long as the 
proper size a/k of a perturbation keeps small in units of 
A kpc, this source can be approximated as 

A galaxy of mass lO^M©, for instance, grew out from 
perturbations of length ^ 10~^ kpc at decoupling; in 
this case the relative driving term on the RHS of ^ is 
initially suppressed by a factor ~ I0~^. This dampening 
is of order ~ I0~^ for galaxies of lO^M©. The full solu- 
tion of (O in our gravitational scenario is plotted in the 
upper panel of Fig. 2 in terms of the fraction Sb/Sd and 
for different masses. 

Using this result, it is straightforward to compute the 
baryon fraction / = (pb/pd) inside shells reaching the 
turn-around. In doing this we use a 7i I power spec- 
trum with F = 0.25, h = 0.7 and erg = 0.84 and we im- 
plicitly assume that linear theory still holds at this epoch 
(at least for its predictions on the ratio Sb/Sd) and that 
the derived / is a good estimate for the baryon fraction 
in the final collapsed object. From the lower panel in Fig. 
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2 one may see that the Yukawa-hke potential acts in such 
a way to avoid that the universal gas fraction {Q.b/^d) 
be locked up in low mass clumps at high redshift. For 
example, we predict / ~ O.SSIb/SI^i in galaxies having 
M = IO^Mq. We note, however, that to fully exploit the 
simplicity of linearized equations, we have neglected the 
statistic of the halo hierarchical merging. 

Our phenomenological assumption lacks an explana- 
tion in terms of fundamental physics. The large absolute 
value of a tends to rule out the interpretation of the 
Yukawa contribution in terms of a scalar field-mediated 
force. In a scalar-field scenario, in fact, the strength pa- 
rameter evaluates a = asaD, where as and an are 
the couplings of the scalar to baryonic and dark matter 
respectively (see e.g. [l^). Since the scalar- baryon cou- 
pling is constrained to be |q!b| < 10~^ by very long base- 
line radio interferometry measurements (see e.g. |27|). a 
dramatic short-range modification of gravity of the order 
of ajj ^ 10'* in the pure dark matter sector could not 
be avoided, at least under the minimal assumption of a 
universal scalar-dark coupling ao- 

In summary, we have argued that the coupling between 
DM and baryons vanishes in the small distance limit and 
it can be effectively described by adding a Yukawa-like 
contribution with universal parameters a — 1,A ^ 1 
kpc to the standard gravitational potential. Our model 
does not affect the relation between dark and visible mat- 
ter on large cosmological scales (i.e. the dynamics of large 
galaxies and clusters) but interestingly it successfully al- 
leviates some inconsistencies between observations and 
CDM theoretical predictions on kpc scales. In particu- 
lar: 

a) the analysis of the rotation curves of LSB galaxies 
shows that the central dark matter density distribution 
predicted under the assumption of a Yukawa potential 
well (/3 — —1.35 ± 0.05) is in excellent agreement with 
the cuspy profile derived in coUisionless CDM simula- 
tions. In particular, the a best fitting values suggest that, 
independently of the adopted Yukawa model, the gravi- 
tational interaction between these two classes of matter 
is totally suppressed in the small distance limit b) this 
potential offers an interesting cosmological mechanism to 
segregate the quantity of baryons that falls into primor- 
dial dark matter perturbations as a function of the fluc- 
tuation mass; this fact may be used to relax estimates of 
the efficiency of various feedback mechanisms proposed 
to suppress the overproduction of dwarf galaxies. 
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